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Abstract 

The presence of nonmagnetic impurities in a 2D "bad metal" depresses the superconduct- 
ing Berezinskii-Kosterhtz-Thouless transition temperature, while leaving the pairing energy 
scale unchanged. Thus the region of the pseudogap non-superconducting phase, where the 
modulus of the order parameter is non-zero but its phase is random, and which opens at 
the pairing temperature is substantially bigger than for the clean system. This supports 
the premise that fluctuations in the phase of the order parameter can in principle describe 
the pseudogap phenomena in high- Tc materials over a rather wide range of temperatures 
and carrier densities. The temperature dependence of the bare superfluid density is also 
discussed. 

PACS: 74.25.-q; 74.40.-hk; 74.62.Dh; 74.72.-h 

1 Introduction 

The differences between the BCS scenario of superconductivity and superconductivity in high- 
Tc materials are well accepted as experimental facts, although there is no theoretical consensus 
about their origin. One of the most convincing manifestations is the pseudogap, or a depletion 
of the single particle spectral weight around the Fermi level (see for example, Q). Another 
transparent manifestation is the temperature and carrier density dependencies of the superfluid 
density in high-Tc superconductors (HTSC) |Q, ^ Q which do not fit the canonical BCS be- 
haviour. In particular the value of the zero temperature superfluid density is substantially less 
than the total density of doped carriers Q. Currently there are many possible explanations for 
the unusual properties of HTSC. One of these is based on the nearly antiferromagnetic Fermi 
liquid model [^j. Another explanation, proposed by Anderson, relies on the separation of spin 
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and charge degrees of freedom. One more approach, which we will follow in this paper, relates 
the observed anomalies to precursor superconducting fluctuations. Different authors argue that 
alternative types of superconducting fluctuations are responsible for the pseudogap, e.g. , and 
the scenario based on the fluctuations of the order parameter phase suggested by Emery and 
Kivelson The latter scenario we believe to be more relevant due to low superfluid density 
and practically 2D character of the conductivity in HTSC mentioned above. A microscopic 2D 
model which elaborates the abovementioned scenario Q has been studied in the papers ||9|, p!o|] . 
The results obtained show that the condensate phase fluctuations indeed lead to features which 
are experimentally observed in HTSC both in the normal and superconducting states Q. It is 
obvious, however, that the present treatment of the phase fluctuations is incomplete due both to 
the oversimplified character of the model and the absence of an explanation for the more recent 
advanced experiments |P, |^, Q on the temperature and doping dependencies of the superfluid 
density. It is well known, however, that the theoretical study of HTSC faces a lot of computa- 
tional difficulties due to, for example, a non-conventional order parameter symmetry, complex 
frequency-momentum dependence of the effective quasiparticle attraction, general form of the 
quasiparticle dispersion law, etc. Therefore, in order to obtain analytical results, we have to 
date only considered nonretarded s-wave pairing in the absence of impurities. (The attempts to 
consider the retardation effects were done in 

Nevertheless, a discussion of the effect of impurities seems to be crucial for a realistic model 
of the HTSC. Indeed, it is known, that the itinerant holes in HTSC are created by doping 
which in turn introduces a considerable disorder into the system, for instance, from the random 
Coulomb fields of chaotically distributed charged impurities (doped ions) [|l2|. Thus one the 
purposes of the present paper is to study the model |lO| but in the presence of nonmagnetic 
impurities. 

In the theory of "common metals" the Fermi energy e p and the mean transport quasiparticle 
time Ttr are independent quantities which are always assumed to satisfy the criterion e pTtr ^ 1 • 
In HTSC which are "bad metals" Q both ep and Ttr are dependent on the doping and the 
abovementioned criterion may fail As an illustration we refer to the remarkable linear 

dependence of the normal state resistivity Q which implies indeed that epTtr may be ~ 1. It 
has been shown [^] that for strongly disordered metallic system superconductivity is absent if 
the scattering-to-pairing ratio exceeds a critical value and the existence of superconductivity in a 
finite range of doping if this ratio is not exceeded. We shall not study this case but rather consider 
here the more usual (and in some sense simple) situation, originally studied in the papers of 



Anderson |13] and Abrikosov-Gor'kov (AG) [14| (see also ||l^), when the superconducting order 
is preexisting and the criterion epTtr ^ 1 is satisfied. 

The Anderson theorem [^] states that in 3D the BCS critical temperature is unchanged in 
the presence of nonmagnetic impurities. However, as discussed in |^, the BCS critical temper- 
ature in 2D is the temperature T* at which the pseudogap opens, while the superconducting 
transition temperature transition is the temperature Tbkt of the Berezinskii-Kosterlitz-Thouless 
(BKT) transition. In contrast to the former, the latter is defined by a bare superfluid density 
(given by the delocalized carriers) which is dependent on (see below) the concentration of im- 
purities. Thus in 2D case the superconducting transition temperature, Tbkt decreases with 

^Of course, the contribution from the phase fluctuations need not be the only or even the major contribution. 
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increasing impurity concentration. 

Thus, in the model under consideration, the relative size of the pseudogap phase, (T* — 
Tbkt)/T* is larger in the presence of impurities than in the clean limit Therefore it can 
be observed over a wider range of densities. The second result obtained is that the value of 
the zero temperature superfluid density is less than the total density of carriers (dopants), so 
that the presence of impurities may contribute into this diminishing and in its turn explain the 
experimental results Q. Finally we attempt to interpret qualitatively the recent experiments 
on the temperature dependence of the superfluid density [Q, ^ within our scenario. 

A brief overview of the paper follows: In Sec. |2| we present the model and derive the main 
equations. In Sec. ^ we compare the results obtained for clean and dirty limits. In particular 
we compare the values of Tbkt , the relative sizes of the pseudogap region and the values of the 
bare superfluid density at T = and T close to Tp. In Sec. ^ an attempt is made to give an 
explanation for the experimental results ^ |^ . 



2 Model and main equations 

Our starting point is a continuum version of the two-dimensional attractive Hubbard model 
defined by the Hamiltonian ^ IC]: 



H 



V2 



+Uimp{r)ipl{x)^a{x) , 



(1) 



where x = r, r denotes the space and imaginary time variables, ipaix) is a fermion field with 
spin a =1, 1, m is the effective fermion mass, /i is the chemical potential, V is an effective local 
attraction constant, and Uimpi^) is the static potential of randomly distributed impurities; we 
take h = ks = 1- The model with the Hamiltonian (|^) is equivalent to a model with an auxiliary 
BCS-like pairing field which is given in terms of the Nambu variables as 



H 



dPrl ^Hx) 



T3 



2m 



T+$(x) - T_$*(x) + T3Uimpir)] ^ (x) + 



|<I>(x) 
V 



(2) 



where t± = (n ± ir2)/2, rs are the Pauli matrices, and ^{x) = V^^x)t-^{x) = 

is the complex order field. Then the partition function can be presented as a functional integral 

over Fermi fields (Nambu spinors) and the auxiliary fields <I>*. 

However, in contrast to the usual method, the modulus-phase parametrization ^(x) = 
p(x) exp[i0(x)] is necessary for the 2D model at finite temperatures (see P, 10, |l^ and refer- 



ences therein). To be consistent with this replacement one should also introduce the spin-charge 
variables for the Nambu spinors 



^{x) = exp[ir3(9(x)/2]T(x) 



(3) 
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with T(x) the field operator for neutral fermions. 

From the Hamiltonian (||), following Q, one can derive an effective one which is the Hamil- 
tonian of the classical XY-model 



where 



HxY = \J{^^,T,p) J d\[ve{x)f 

T °° f 

n=—oo 



(4) 



(5) 



is the bare (i.e. unrenormalized by the phase fluctuations, but including pair breaking thermal 
fluctuations) superfluid stiffness. Here 



with 



(^(ia;„,k)) 



K + +e2(k) 



(6) 



1 + 



7r(2n + l)r 



(7) 



is the AG |15] Green's function of neutral fermions averaged over a random distribution of 
impurities and written in the Nambu representation |17, |l8|. In writing ^ we assumed that 
k)^(iu;„, k)) — (^/(^'^n, k)) (^(ia;„, k)). This approximation, as shown by AG ||l^, does 
not change the final result for J. Note also that the Green's function (|^) is valid only when 
GpTtr ^ 1 which demands the presence of a well developed Fermi surface, which in turn implies 
that ;U ~ e_p. Thus one cannot use the expression @ in the so called Bose limit with |U < [^. 
On the other hand, Fermi surface can be formed even in the bad metals when loffe-Regel-Mott 
criterion proves to be fulfilled [p^]. 

Substituting @ into (|5|) , and using the inequalities /i S> T, /) to extend the limits of integra- 
tion to infinity, one arrives at 



An An ^ 



dx 



(8) 



Eq.(|8|) is formally divergent and demands special care due to the fact that one has to perform 
the integration over x before the summation |15|. Finally one can formally cancel the divergence 
[15 1 to obtain 



J 



E 



(9) 
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The temperature of the BKT transition for the XY-model Hamiltonian is determined by 
the equation 

TT 

Tbkt = ■^-^(/",7bkt,p(^,7bkt))- (10) 

The self-consistent calculation of Tbkt as a function of the carrier density Uf = me p /tt requires 
additional equations for p and p, which together with form a complete set 

When the modulus of the order parameter p(x) is treated in the mean field approximation, 
the equation for p takes the following form 

2p _ ^ f dPk 



I (^tr[ri(g(it^„,k))], (11) 
which formally coincides with the gap equation of the BCS theory. This coincidence allows 



V ^ 



one to use the Anderson theorem |13| which states the dependence of /o(r) is the same as that 
for the clean superconductor and not affected by the presence of nonmagnetic impurities. It is 
important to recall that this theorem is, of course, valid only for the s-wave pairing and small 
disorder. 

There is, however, both physical and mathematical differences between the gap in the BCS 
theory and p |^ 10 1. In particular the temperature Tp which is estimated by the condition /j = 
is not related to the temperature of the superconducting transition, but is interpreted as the 
temperature of pseudogap opening T* (see details in [^). The main point, which we would like 
only to stress here, is that due to the Anderson theorem ||T^ the value of Tp does not depend 
on the presence of impurities, while the temperature Tbkt, as we will show, is lowered. 

The chemical potential p is defined by the number equation 

/(02t^N(a(iu;n,k))] =nj. (12) 
Since we are interested in the high carrier density region the solution of (O) is // ~ e^?, so that 



— oo ' 



in Eqs.dlD - (|llD one can replace p by tp. 

Having the temperatures Tp and Tbkt as functions of the carrier density one can build the 
phase diagram of the model which consists of three regions. The first one is the supercon- 
ducting (here BKT) phase with /o 7^ at T < Tbkt- In this region there is algebraic order, or 
a power law decay of the correlations. The second region corresponds to the pseudogap 

phase (Tbkt < T < Tp). In this phase p is still non-zero but the correlations mentioned above 
decay exponentially. The third is the normal (Fermi-liquid) phase at T > Tp where p = 0. Note 
that (<^(x)) = everywhere. While the given phase diagram was derived for the idealized 2D 
model, there are indications that even for as complicated layered systems as HTSC the value 
of the critical temperature for them may be well estimated using Tbkt i^] even though 
the transition undoubtedly belongs to the 3D XY class. It was also pointed out in |jl9[ that a 
nonzero gap in the one-particle excitation spectrum can persists even without long-range order. 
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3 The comparison of clean and dirty limits 



3.1 Clean limit 

The transport time Ttr is infinite in the clean hmit, so that 



Ji€F,T,p{eF,T)) 



epp'T 



(u;2+p2)3/2- 



(13) 



Near Tp one can obtain from ( |13| 



J{eF,T ^ Tp ,p ^0) = -TTTJTT^^F, 



16^3 



(14) 



where C{^) is the zeta function. This expression must coincide with the result from Q which was 
derived using the opposite order for the summation and integration. Inserting the well-known 
dependence of p{T) (see, for example [pi] ]) 



p\T^T-) 



'-'-T^d-^ 



7C(3)' 



(15) 



and then substituting (jl^) into (^0|) one obtains the following asymptotic expression for the 
BKT temperature in the clean limit for high carrier densities [22,^, 23 1 



BKT 



1-^1, Tbkt < 



In the high density limit one can also use the equation 



7 



ep, 



(16) 



(17) 



where 7 ~ 1.781 and is the energy of the two-particle bound state in vacuum which is a more 



convenient parameter than the four-fermion constant V |24, y. 

It is obvious from (^) and (^^ that the pseudogap region shrinks rapidly for high carrier 
densities and one may ask (see, for example, |25]) whether this scenario can explain the 
pseudogap anomalies which are observed over a wide range of temperatures and carrier densities, 
since in the clean limit the relative size of the pseudogap region {Tp — Tbkt)/?)? is, for instance, 
less than 1/2 when the dimensionless ratio ep/\ei,\ < 1287^/7r^ ~ 41. A crude estimate for the 
dimensionless ratio for optimally doped cuprates gives ep/\eh I ~ 3 • 10^ - 10^ [|6| which indicates 
that in the clean superconductor the pseudogap region produced by the phase fluctuations is 
too small. Of course, all these estimations are qualitative due to the simplicity of the model. 

The value of the bare superfluid density, ns{T) is straightforwardly expressed via the bare 
phase stiffness, ns{T) = AniJiT). In particular, it follows from (13) that ns{T = 0) = Uf. This 
is not surprising since ns{T = 0) must be equal to the full density Uf for any superfluid ground 
state in a translationary invariant system |27] and the clean system is translationary invariant. 



In 2D for s-wave pairing the high density limit is in fact equivalent the weak coupling BCS limit. 
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We note, however, as stated above that in HTSC UgiT = 0) nf Q. Substituting (^) into 
(13) one obtains for T close to Tp the bare superfluid density as ns{T — > T~) = 2nj(l — T /Tp). 
This behaviour of the bare superfluid density is formally the same as the behaviour of the full 
superfluid density in the BCS theory. Nevertheless it is important to remember that the full 
superfluid density in the present model undergoes the Nelson-Kosterlitz jump at Tbkt and is zero 
for T > Tbkt- We note that one can probe experimentally both the bare superfluid density in 
high-frequency measurements [0] and the full superfluid density in low-frequency measurements 



3.2 Dirty limit 

In the dirty limit the quasiparticle transport time Ttr is small {rtr *C p^^iT = 0)) so that one 
can neglect the radical in the bracket of @ ||l^ . The remaining series is easily summed and one 
obtains for the bare superfluid stiffness 

J{eF,T,p{eF,T),Ttr) = ^^^''^ tanh^. (18) 



As explained above, due to the Anderson theorem, the expressions (|T5|) for p and ( p^ ) for Tp 
remain unchanged in the presence of impurities. Again substituting (|15|) into ( p!8|) one obtains 



tBKT 



Tp ( 1 ) , Tbkt < Tp. 



(19) 



One can see that the size of the pseudogap region is now controlled by the new phenomenological 
parameter Ttr which is an unknown function of ep for HTSC. The experimental data [||] suggest 
that Ttr is almost independent on doping level in the underdoped region. 

It is difficult to obtain more than a qualitative estimate using Eq.([T9|) since in its derivation 
we have assumed that epTtr ^1. In HTSC however, as discussed above (see also [|l2|), this 
assumption is not always justified. Bearing in mind that the dirty limit implies that the condition 
T^^ ^ p{T = 0) ^ Tp is satisfied, one can easily see that the value of Tbkt for this case is less 
than that given by (16) for the clean superconductor. Since impurities are inevitably present in 



HTSC, phase fluctuations can in fact give rise to a pseudogap region that is of comparable size 
to that experimentally observed. We note that our arguments are in fact quite similar to that 
given in [p^ j for the best observing conditions for the BKT physics in superconducting films. 
However, in contrast to this paper, the gap opening below Tp is particularly emphasized here. 
While Eq.(19) was derived under assumption Tbkt ^ Tp in general case when Tbkt can be 



substantially less than Tp the self-consistent equation (px|) with J(ei?, Tbkt, /o(eF) Tbkt)) given 



by (18) must be solved. Recall, however, that to make any quantitative estimates, the more 



realistic d-wave model has to be considered and the inequality epTtr J?> 1 should not be assumed 
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The value of the zero temperature superfluid density is now given by ns{T = 0) = irnfTtrP <C 
Uf since TtrP <C 1. This does not contradict the results of |27] because the system is not 
translationary invariant in the presence of impurities |28|. Furthermore as one can see, the 
low value of the superfluid density in HTSC may be related to the impurities which are 
inevitably present in HTSC. Another reason which leads to lowering of the superfluid density is 
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the presence of lattice which also destroys a continuous translational invariance. We note that 
as was pointed [^] quantum fluctuations also lead to a decrease in the superfluid density. 



4 The temperature dependence of the bare superfluid density 

In this section we try to correlate the temperature dependence of the observed in-plane resistivity 
Pab(T) with the recently measured temperature dependence of the bare superfluid density |||]. 

For T > Tbkt the expression for the bare superfluid density in the dirty limit (|^) can be 
rewritten in terms of the in-plane conductivity, a = e^nfTtr/m, where e is the charge of electron: 

Jiai€F,T),p{eF,T)) = ^^tanh^. (20) 

The in-plane resistivity, pab ~ in cuprates has been extensively studied and its 
temperature and concentration dependencies must reflect the pseudogap properties observed in 
other experiments. One can say that Pab{T) is linear above T* ~ Tp and roughly linear between 
Tbkt and Tp but with a lower slope. Thus in the interval Tbkt < T < Tp the resistance can 
be approximately written as Pab{T) = aT + b, where a and b are functions of ep but not of 
temperature. 

Now, substituting a ~ P^h^T) into Eq.(pO[), one obtains 

UsiT) ~ — ^ tanh (21) 
^ ' aT + b 2T ^ ' 

Our estimations based on Eq.(21) are shown in Fig. ||. One can see that, in contrast to the 
almost linear BCS dependence of ns(T), we have convex behaviour and the superfluid density 
becomes zero at Tp. We stress that the curvature of ns(T) is the result of both the temperature 
dependence of p{T) and o"(T) for Tbkt < T < Tp. More importantly the slope of the curve 
ns(T) at Tp for the dirty metal is substantially less than for the clean one. The experiment |Q 
shows the same curvature for ns{T) but indicates that the bare superfluid density disappears 
at a lower temperature, Tg < T* . Since the slope dns{T)/dT at Tp is very small, as predicted 
by Eq.(pl|) and observed experimentally, the non-zero value of ns{T) between T^ and Tp may 
however simply be too small to be experimentally observed. A definitive answer to this question 
demands further experiments and theoretical studies. In particular, ^-fluctuations should be 
taken into account |^, |l^ . 

One can also comment on the experimentally observed change in the curvature of the full 
superfluid density, J\fs{T), with changing carrier density even though A/'s(T) cannot be directly 
linked to the bare superfluid density n<j(T) discussed here. Although the full superfluid density 
disappears above Tbkt, the curvature present in the bare superfluid density ns{T) seems to be 
retained as a curvature in the full superfluid density, 7Vs(T), below Tbkt ||2|, ^. For low carrier 
densities (the underdoped region) the pseudogap region, Tbkt < T < Tp, is larger and therefore 
the curvature in ns{T) is more pronounced. This behaviour seems to be reflected in the full 
superfluid density, 7Vs(T) below Tbkt Q. It is important however to study experimentally and 
theoretically the concentration dependence of the bare superfluid density, ns(T) to make a full 
comparison with the results from for (T) . 

The experimental data of S also show that Ms{T) does not display the Nelson-Kosterlitz 



jump. This is probably related to the influence of the interlayer coupling, see Refs. in |10|. 
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5 Conclusion 



Since in HTSC the pairing scale T* is different from the superconducting transition temperature 
the role of non-magnetic impurities is not traditional and they in fact define the superconducting 
properties of a "bad metal". In particular the presence of nonmagnetic impurities strongly 
increases the size of the pseudogap phase originating from the fluctuations of the phase of the 
order parameter. In addition the behaviour of the superfluid density in the presence of the 
impurities is closer to that experimentally observed. 

Our results are only qualitative since we have considered the model with non-retarded s- 
wave attraction and an isotropic fermion spectrum. However, it is likely that for d-wave pairing, 
the properties obtained will persist. There is, of course, a problem why a strong disorder does 
not destroy d-wave superconductivity when non-magnetic impurities are pair-breaking. As was 
suggested by Sadovskii |Q even d-wave pairing may persist if coupling is strong enough. Further 
studies are necessary, for example, it is important to explain the concentration dependence of 
the superfluid slope, dns{T) / dT at T = ^. Our results also indicate that it would be 
interesting to study the BCS-Bose crossover problem in the presence of impurities, especially in 
d-wave case [pO| . 

We gratefully acknowledge V.P. Gusynin and Yu.G. Pogorelov for many useful discussions. One 
of us (S.G.Sh) is grateful to the members of the Department of Physics of the University of 
Pretoria for hospitality. R.M.Q and S.G.Sh acknowledge the financial support of the Foundation 
for Research Development, Pretoria. 
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